Evidence of Luttinger liquid behavior in one-dimensional dipolar quantum gases 
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The ground state and structure of a one-dimensional Bose gas with dipolar repulsions is investi- 
gated at zero temperature by a combined Reptation Quantum Monte Carlo (RQMC) and bosoniza- 
tion approach. A non trivial Luttinger-liquid behavior emerges in a wide range of intermediate 
densities, evolving into a Tonks- Girardeau gas at low density and into a classical quasi-ordered state 
at high density. The density dependence of the Luttinger exponent is extracted from the numerical 
data, providing analytical predictions for observable quantities, such as the structure factor and 
the momentum distribution. We discuss the accessibility of such predictions in current experiments 
with ultracold atomic and molecular gases. 
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The realization of Bose-Einstein condensation (BEC) 
in trapped ultracold quantum atomic gases [l[ is at the 
frontier of modern atomic and molecular, optical and 
condensed- matter physics Q. Especially fascinating ex- 
perimental results arise from the possibility of tuning the 
atomic interactions. Use of Fano-Feshbach resonances Q 
to change magnitude and sign of the s-wave scattering 
length a characterizing the contact interactions, has al- 
lowed e.g. the observation of collapsing Bose conden- 
sates [B| and of the crossover from a BEC to a Bardeen- 
Cooper-Schrieffer-like transition Q . 

More recent experiments have demonstrated that the 
range of the interactions can also be manipulated. Dipolc 
interactions with long-range anisotropic character have 
been observed in 52 Cr atoms [7j after exploiting the large 
magnetic moments of this atomic species, that is [id rs 
6/ib with fiB being the Bohr magneton. A BEC contain- 
ing up to 50 52 Cr atoms has then been obtained below 
a transition temperature T c ~ 700nK |8| and its dynam- 
ical behavior is being investigated [9(. Promising pro- 
posals to tune and shape the dipolar interaction strength 
in quantum gasees of heteronuclear polar molecules have 
more recently been suggested [l(| ■ Significant theoretical 
predictions have accompanied such realizations [Til ] . The 
stability diagram of anisotropic confined dipolar gases 
has been predicted to be governed by the trapping ge- 
ometry 
studies 



known to emerge from the collective character of the 
single-particle degrees of 

freedom, despite the absence of broken symme- 
tries 21] . Bosons are known to arrange in a Luttinger- 



12 , [13[ , as corroborated by Path-Integral QMC 
14]. Different conclusions are reached by more 



recent Diffusion QMC simulations including the depen- 
dence of a on the dipole interaction [HI . 

Tuning of the interactions can be combined with 
the enhancement of quantum fluctuations after reducing 
their dimensionality by e.g. storing them in elongated 
traps [l(| [13], which could be relevant to applications 
such as precision measurements [l8j], quantum comput- 
ing [Hi ] , atomtronic quantum devices, and theoretical in- 
vestigations of novel quantum phase transitions [13] ■ 

In the case of quasi one-dimensional (ID) condensates 
with short-range interactions, a rich phenomenology is 



liquid state, with single particles being replaced by col- 
lective density excitations [H, [23| . Strong repulsion may 
also lead to the fermionization of interacting bosons in 



the so-called Tonks-Girardeau (TG) regime [24 l25l. |2q|. 
Experiments in elongated traps have provided evidence 
for such ID fluctuations [16j . 

In the case of quasi-lD condensates with dipolar inter- 
actions, an interesting question arises whether the quan- 
tum fluctuations are sufficiently enhanced to drive the 
BEC in a strong- coup ling regime. More recent Diffusion 
QMC simulations [27J for a homogeneous ID dipolar Bose 
gas have revealed a crossover behavior with increasing 
linear density, from a liquid-like superfluid state to an 
ordered, normal, state with particles localized at lattice 
sites. For low values of the density, the system is seen to 
behave as a Tonks-Girardeau gas [241 ] . 

A consistent theoretical description of the whole 
crossover region would be useful to predict the behav- 
ior of observable quantities in the experiment. It is 
well known [22j] that the low-energy properties of a one- 
dimensional system of interacting bosons can be ob- 
tained from an effective quadratic Hamiltonian within a 
Luttinger-liquid theory, provided that the interaction is 
sufficiently short ranged. Besides, ID fermionic models 
with interactions falling off as Jjx a have been studied 
in the past by many authors 35, [3^], concluding that 
for a > 1 the long-range part of the interaction decays 
rapidly enough to preserve the linear energy dispersion 
ui(k) = u\k\ typical of short-range models. 

In this Letter we provide evidence for a robust 
Luttinger-liquid behavior of a ID dipolar Bose gas in 
the crossover region, evolving into a Tonks-Girardeau gas 
at low density and into a classical quasi-ordered state 
at high-density. We reach this conclusion by comparing 
theoretical results from a bosonization approac h ag ainst 
Reptation Quantum Monte Carlo simulations [28(. By 
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analyzing the scaling behavior of the static structure 
factor expected for a Luttinger liquid, we find small 
Luttinger exponents at intermediate-to-high density val- 
ues, signalling strong repulsions at short distances. The 
knowledge of Luttinger parameters is then exploited to 
determine observable properties of the fluid in a param- 
eter range accessible to experiments. To the best of our 
knowledge this is the first demonstration of such strongly 
correlated liquid behavior in atomic Bose gases with dipo- 
lar interactions [29I ]. 

We consider N atoms with mass M and permanent 
dipoles moments arranged along a line. For the purposes 
of the present work, we assume a — > after e.g. ex- 
ploiting a Fano-Feshbach resonance. In effective Rydberg 
units Ry* = h 2 j(2Mr\) the Hamiltonian is 



^ ( T 2 ^ ; r ^ ^ / 



1 



(1) 



where we have defined the effective Bohr radius ro = 
MCdd/(2Trh 2 ) in terms of the interaction strength Cdd- 
One has Cdd — MoMd anc ^ ^dd — d 2 /eo for magnetic 
and electric dipoles respectively, where Hd and d are the 
magnetic and electric dipole moments and a and eo are 
the vacuum permittivities. The dimcnsionless parameter 
governing the system is r s = l/(nro), with n being the 
linear particle density and r s ro the unit length. 

We analyse the low-energy structural properties of the 
fluid. They can be accessed in experiments with atomic 
gases by means of e.g. Bragg scattering techniques. Here, 
we infer them from QMC numerical simulations of the 
static structure factor S(k) on a box with size L and 
periodic boundary conditions. In terms of the density 
operator n this is 



S(k) = 



dxe 



— ika 



(n(z)n(O)). 



(2) 



An accurate size-scaling analysis of S(k) is needed to 
our purposes. To this aim we resort to the Reptation 
Quantum Monte Carlo of Baroni and Moroni [28| , which 
in essence is a path-integral technique at zero temper- 
ature giving direct access to the ground-state wavefunc- 
tion. This makes the evaluation of the ground-state prop- 
erties conceptually simple and the extraction of the corre- 
lation functions in imaginary time immediate [3^|, with 
the additional possibility of determining the excitation 
spectrum. We use a trial wave-function made up of a two- 
body Jastrow factor iptriai(R) = ^h<j exp[u(\x{ — Xj\)]. 
Removal of the divergences in Hiptriai/ 'Atrial due to the 
potential energy in the x — *■ limit, is obtained after 
imposing the cusp condition, which yields u(x — > 0) — 
-^[ufx~ & 

The RQMC data for S(k) are reported in Fig. Q] for 
N = 40 and different values of nr = 0.01,50,100 
and 1000. Free fermion-like behavior, typical of the 




k/n 

FIG. 1: S(k) in dimensionless units for a dipolar gas with 
N = 40 particles and different values of nro = 0.01, 50, 100 
and 1000. Decreasing slopes as k — > and the emergence of 
additional peaks correspond increasing nro values. 



Tonks- Girardeau regime, is seen at low density values 
(nro = 0.01 in the Fig.[T]). Quasi-Bragg peaks emerge at 
reciprocal lattice vectors k/n — 2nm (m integer). Their 
number increases with increasing nro, evidencing the oc- 
currence of a quasi-ordered state. Investigation of the 
RQMC data at different values of N show that the peak 
heights scale as N Um with a m < 1 being an exponent 
depending on to, while their intensity diverges with the 
system size. 

Both these characteristics are reminiscent of Luttinger- 
liquid behavior and therefore we turn on analysing the 
data within a bosonization description. In terms of the 
conjugate bosonic operators 4> and n, the Hamiltonian 



for a Luttinger liquid reads [21 1 



H 



uKinTl) 2 + -{d x 

A 



(3) 



where K is the Luttinger exponent and u the sound veloc- 
ity. Knowledge of K and u allows analytical expressions 
for physical quantities of interest. The structure factor 
S(k) can be analytically calculated inserting in Eq. @ 
the bosonized expression of the density operator n. that 
is h(x) = -^d x <t> + Y, m ^o A mne 2i ( m ^-^ 0, HI, 
where A m are non-universal dimcnsionless constants. In 
dimensionless units for k, S(k) turns out to be So(k) = 
Kkj (27r) for k — > 0. For k close to k = 27tto on the scale 
of the cutoff parameter Tr/a, we find [37| : 



S(k) 



J2 A 2 m n 2 ^ m (k±2^m) , 

m^O 



(4) 



where H m (k) is explicitly given as 

£ m (fc) - L(l - e- 2rra l L ) 2m2K a(k- m 2 K, L) (5) 
x 2 F 1 (m 2 K, m 2 K + \k\L/(2n); 1 + |fc|L/(27r); 



where cr(fc; to 2 = r[TO 2 A'+|/c| J L/(27r)][r(TO 2 A')r(l- 
fc|L/(27r))] _1 , 2F1 the hypergeometric function. 
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The low-fc linear behavior of S(k) can be fitted to 
extract the Luttinger K at different densities. The re- 
sulting K(n) is displayed in Fig. O The inset shows 
S(k — » 0) for selected densities, where descreasing slopes 
correspond to increasing values of nr = 50, 100, 1000, 
the symbols representing the RQMC data and the lines 
the Luttinger prediction S(k — > 0) = Kk/(2n). K re- 
mains below the Tonks-Girardeau limit in a wide range 
of densities, signalling the occurrence of strong correla- 
tions. As nro < 1, K tends to the Tonks-Girardeau 
limit K = 1. In the opposite high-density limit, K 
tends to vanish following the classical-theory prediction 
K(n -> oo) ~ 7r(6C(3)nr )- 1 / 2 . 
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FIG. 2: Luttinger exponent K obtained from fitting the 
RQMC data for S(k) in the k — > limit. For comparison, 
the curves K(n — > 0) = 1 and K(n —* oo) = 7r(6C(3)nro)~ 1/ ' 2 
are also displayed as dashed lines, corresponding to the Tonks- 
Girardeau nro — * limit and to the classical prediction for 
nro — > oo. The solid line is a guide to the eye. The inset 
zooms on the low-fc portion of S(k) at selected nro (see text). 

S m (fc) acquires a simplified form depending on 
whether 7 = m 2 K < 1/2 or 7 > 1/2 as L — > 00. For 
7 < 1/2 we get E m (fc) - L 1 - 2m " K sin{nm 2 K)f3(K] L) 
with (3{K- L) = T(l - 2m 2 K)T(m 2 K + kL/(2w))/T(l ~ 
m 2 K + kL/(2n)), indeed revealing the presence in S(k) 
of quasi-Bragg peaks at |fc| = 27rm and a scaling behav- 
ior in the vicinity of the peak. The peak intensities di- 
verge with the system size L and their number increases 
with decreasing K (namely increasing n). The peak's 



height scales as N 



l-2m 2 K 



For 7 > 1/2 instead we get 



E m (fc 0) ~ 2 1 - 2 T[r(2 7 - l)][r( 7 ) 2 ]- 1 and the peaks 
disappear. 

Besides the linear k behavior, evidence of a Luttinger 
liquid comes also from the scaling behavior ofS(k). Fig. [3] 
displays the comparison between the RQMC S(k) close 
to k = 2n towards the quasi-ordered region with nro — 
50, and the predictions of the Luttinger-liquid theory 
Eqs. Q-© with m = 1. We actually plot S(k)/S(2nn) 
for different N values and verify that they indeed collapse 



on a single curve. This signature of Luttinger-liquid be- 
havior against the RQMC data is robust in the whole 
crossover region. This is the central result of the present 
work. Similar results are indeed found at the remaining 
density values above titq = 1, while the scaling is absent 
below nro = L The same scaling analysis holds for the 
second peak at k — Air when it is present. At k = 67r 
(m = 3), we have 2m 2 K > 1 for K < 0.055 and the 
divergent peak is seen to show up in the simulation only 
for the data set with nro = 1000, as expected. 

The peak's asymmetry visible in Fig. [3] results from 
band curvature effects. Asymmetric structure factors 
are for example obtained within exact calculations on 
Calogero-Sutherland models [32} which indeed possess a 
non-linear spectrum, and could be also viewed as anhar- 
monic phonon effects 
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FIG. 3: S(k)/S(2ir) in n" 1 units vs. iV(fc-27r)/40 at nr = 50 
for different N values (see legend). Symbols represent RQMC 
data, the solid line is the Luttinger-liquid (LL) prediction. 

As a consistency check, we have verified that the 
K values obtained from the RQMC S(k) agree with 
those determined from the RQMC energy per particle 
e, as follows. Since Hamiltonian ([1]) is invariant under 
Galilean boostsfUl, we have uK = him/M. Within 



the bosonization procedure [21|, [22j , the compressibility 
is x — KChnun 2 )" 1 . Comparing this result with the def- 
inition of x m terms of density derivatives of e, that is 
X^ 1 = n 2 (ne(n))" ', we obtain K = tvK^J n/[m(ne(n))"]. 
This result stems from the trivial fact that the exact 
RQMC satisfies the compressibility sum-rule. However, 
we have also demonstrated on a purely formal basis that 
our Luttinger-liquid scheme satisfies the x sum-rule too. 

The knowledge of K{n) allows us to determine other 
observables such as the momentum distribution n(fc), 
which can be accessed in current atomic BEC's experi- 
ments via e.g. the analysis of time-of-flight images. Since 
in the continuum the boson creation operator is repre- 
sented by ipg{x) ~ exp[— iir J x dyll(y)} [2lj], we find that 
n(k) oc k x l 2K ~ x as L — » 00. As the density increases, the 
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divergence of n(k) at k = is expected to be gradually 
reduced until it disappears for K < 1/2. 

In conclusion, we have shown that dipolar gases in a re- 
duced ID geometry arc promising candidates to observe 
a strong-coupling Luttinger- liquid behavior with K < 1. 
The Luttinger-liquid behavior is robust over a wide range 
of nrg values, evolving into a Tonks-Girardeau gas at low 
density and into a classical quasi-ordered state at high- 
density. Further testing of the Luttinger-liquid behav- 
ior can be provided by the excitation spectrum in the 
crossover [391 ] where, as from preliminary data, the same 
combination of bosonization and RQMC techniques can 
be successfully applied. Observation of the Luttinger liq- 
uid in the crossover region is within the reach of current 
and future experiments in elongated traps. The effec- 
tive ID character of trapped dipolar gases is governed 
by the condition rT 1 3> l±, with l± = (ro/(4a\)) 1 ' 5 a± 
ensuring the suppression of forward scattering [10| . Val- 
ues of the transverse oscillator length a± — yJh/(Mu>±) 
related to the harmonic angular frequency lo±_ can be 
pushed to 50 nm. Then, dipolar SrO molecular gases in 
elongated traps are promising candidates, as we might 
have l± ~ 0.2/i m while ro w 240 /im, yielding access 
to nro values up to the quasi-ordered regime. For dipo- 
lar 52 Cr gases instead we can have l± rs 31 nm while 
ro ~ 4.8 nm, yielding access to nr$ values remaining 
down in the Tonks-Girardeau regime. If ultracold SrO 
dipolar quantum gases were to be created in such quasi- 
1D conditions, we predict the observation of a power-law 
behavior of the peak heights in the structure factor and a 
flattening of the momentum distribution with increasing 
nrg as fingerprints of Luttinger-liquid behavior. 

We thank Dr. Paolo Pedri for very useful suggestions 
on the quasi-lD realizations. RC, SDP and MLC thank 
Prof. G. C. La Rocca for useful discussions. SDP is 
indebted to Dr. S. Moroni for helpful discussions and 
QMC support. 



[1] 



[2] 



[3] 

[4] 
[5] 

[6] 



[7] 



M. J. Anderson et al, Science 269, 198 (1995); K. B. 
Davis et al, Phys. Rev. Lett. 75, 3969 (1995); C. C. 
Bradley et al., ibid. 75, 1687 (1995). 
Bose-Einstein Condensation in Atomic Gases, Proc. 
Intl. School of Physics "Enrico Fermi", Course CXL, M. 
Inguscio, S. Stringari, and C. Wieman Eds. (IOS Press, 
Amsterdam, 1999). 

H. Feshbach, Ann. Phys. 5 357 (1958); U. Fano, Phys. 

Rev. 124 1866 (1961). 

S. Inouye et al, Nature 392, 151 (1998). 

J. Roberts et al, Phys. Rev. Lett. 86, 4211 (2001); T. 

Weber et al, Science 299, 232 (2003). 

C. A. Regal et al, Phys. Rev. Lett. 92, 040403 (2004); 

M. Bartestein et al, Phys. Rev. Lett. 92 120401 (2004); 

M.W. Zwierlein et al, Phys. Rev. Lett. 92, 120403 

(2004); T. Bourdel et al, Phys. Rev. Lett. 93 050401 

(2004); M.W. Zwierlein et al, Nature 435, 1047 (2005). 

J. Stuhler et al, Rev. Lett. 95, 150406 (2005). 



[8 

[9 
[10 

[11 
[12; 
[13 



[14 
[15 

[16 



[17 

[is; 

[19 

[20 

[21 

[22 
[23 

[24 
[25 
[26 

[27 

[28 

[29 

[30 
[31 

[32 
[33 
[34 

[35 

[36 
[37 

[38 
[39 



A. Griesmaier et al, Phys. Rev. Lett. 94, 160401 (2005). 

S. Giovanazzi et al, Phys. Rev. A 74, 13621 (2006). 

H. P. Buchler et al. , |cond-mat /0607294| (2006) . 

See e.g. M. Baranov et al, Phys. Scripta T 102, 74 

(2002). 

L. Santos et al, Phys. Rev. Lett. 85 1791 (2000); K. 
Goral et al, Phys. Rev. A 61, 05160RR) (2000). 
The use of time-dependent external fields for stability has 
been suggested by S. Giovanazzi et al, Phys. Rev. Lett. 
89, 130401 (2002). 

K. Nho et al, Phys. Rev. A 72, 023615 (2005). 
D.C.E. Bortolotti et al, jcond-mat/060 4432 (2006). 

D. Hellweg et. al., Acta Phys. Pol. B, 73, 781 (2001); 
A.Gorlitz et al, Phys. Rev. Lett. 87, 130402 (2001); S. 
Richard et al, Phys. Rev. Lett., 91, 010405 (2003). 
For a collection of papers: Quantum Gases in low di- 
mensions, Pricoupenko et al. Eds., J. Physique IV 116 
(2004). 

See e.g. MA. Kasevich, Science 298 1363 (2002); I. Caru- 
sotto et al, Phys. Rev. Lett. 95 093202 (2005). 
See e.g. I. Bloch, Nature Physics 1 23 (2005) and refer- 
ences therein. 

K. Goral et al, Phys. Rev. Lett. 88, 170406 (2002); H. 
Pu et al, Phys. Rev. Lett. 87, 140405 (2001); H. Pu et 
al. |cond-mat/0203066| (2003) . 

See e.g. T. Giamarchi, in Quantum Physics in one Di- 
mension, (Oxford University Press, Oxford, UK, 2004). 

F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

H. Moritz et al, Phys. Rev. Lett. 91, 250402 (2003); 

M.A.Cazalilla, J. Phys. B: At. Mol. Opt. Phys. 37, 

Sl(2003); D.Petrov et al. J. de Phys. IV 116, 3 (2004). 

M. Girardeau, J. Math. Phys. 1, 516 (1960). 

T. D.Schultz, J. Math. Phys. 4, 666 (1963). 

On the possible realization of the TG regime see B. Pare- 

des et al, Nature 429, 277 (2004). 

A. S. Arkhipov, G. E. Astrakharchik, A. V. Belikov, and 

Yu. E. Lozovik, JETP Lett. 82, 39 (2005). 

S. Baroni and S. Moroni, Phys. Rev. Lett. 82, 4745 

(1999). 

Similar results are known in lattice models with hard- 
core [3(| or continuum models with 1/r 2 repulsions [3ll . 

m. 

G. Gomez-Santos, Phys. Rev. Lett. 70, 3780 (1993); ibid 
73, 1053(E) (1994). 

F. Calogero, J. Math. Phys. 10, 2191 (1969); ibid 2197 
(1969); B. Sutherland, J. Math. Phys. 12, 246 (1971); 
ibid 251 (1971). 

E. R. Mucciolo, B. S. Shastry, B. D. Simons, and B. L. 
Altshuler, Phys. Rev. B49, 15197 (1994). 

S. De Palo, S. Conti and S. Moroni, Phys. Rev. B 69, 
035109 (2004). 

We obtain the same results also using either the Luttinger 
form u(x) oc sin a (7r:r/L) or gaussians centered on the 
lattice sites as trial wavefunctions. 

Y. Tsukamoto and N. Kawakami, J. Phys. Soc. Jpn. 69, 
149 (1999). 

H. Inoue and K. Nomura, J. Phys. A 39, 2161 (2006). 
For further details we refer to E. Orignac et al, in prepa- 
ration. 

H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 

S. Giovanazzi et al, cond- mat /041 1678 (2004) predict a 

roton spectrum within a quasi-lD mean-field calculation. 



